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We analyze the low energy properties of a device with N + 1 quantum dots in a star configuration. 

A central quantum dot is tunnel coupled to source and drain electrodes and to N quantum dots. 
Extending previous results for the N = 2 case we show that, in the appropriate parameter regime, the 
low energy Hamiltonian of the system is a ferromagnetic Kondo model for a S = (N — l)/2 impurity 
spin. For small enough interdot tunnel coupling, however, a two-stage Kondo effect takes place as 
the temperature is decreased. The spin 1/2 in the central quantum dot is Kondo screened first and 
at lower temperatures the antiferromagnetic coupling to the side-coupled quantum dots leads to an 
underscreened S = N/2 Kondo effect. We present numerical results for the thermodynamic and 
spectral properties of the system which show a singular behavior at low temperatures and allow to 
characterize the different strongly correlated regimes of the device. 


I. INTRODUCTION 

The ferromagnetic and underscreened Kondo models 
are examples of quantum impurity problems where the 
low energy properties cannot be described using Fermi- 
liquid theory mu. In these models, a local magnetic 
moment decouples asymptotically at low energies from 
a non-interacting electron reservoir leading to a singular 
behavior in the thermodynamic, dynamic, and transport 
properties of the impurity Pi. The resulting asymptot¬ 
ically free magnetic moment is extremely sensitive to a 
Zeeman splitting and at zero temperature is easily polar¬ 
ized by any non-zero external magnetic field. The lifting 
of the degeneracy by the magnetic field changes the na¬ 
ture of the ground state driving the system into a Fermi 
liquid regime. This extreme sensitivity to external fields 
could be used in nanoscopic devices to control the elec¬ 
tronic and thermal transport and to generate spin cur¬ 
rents [0 • 

Underscreened Kondo physics has been recently re¬ 
ported in spin-1 molecular junctions and multilevel quan¬ 
tum dots where a single channel (not counting spin) in 
the electrodes is relevant for the screening process for a 
wide range of temperatures jTTHl4] . In these systems, 
however, the magnetic moment is expected to be fully 
screened at low enough temperatures since the coupling 
to a second conduction channel cannot be ruled out by 
symmetry considerations. Quantum dot (QD) devices 
in semiconductor heterostructures offer the possibility of 
tuning the parameters and setting the geometry to obtain 
different quantum impurity models m- Devices having 
three QDs in a star configuration have been shown to 
lead to the spin-1/2 ferromagnetic Kondo model TSH1 Oj . 
The realization of such devices would allow the first un¬ 
ambiguous measurement of ferromagnetic Kondo physics. 
Here we show that including additional side-coupled QDs 
the system can be used to construct, by appropiatedly 
setting the parameters, the ferromagnetic Kondo model 
for a spin S = (N — l)/2, where N is the number of side- 
coupled QDs, and the underscreened Kondo model for a 
spin S = N/2. 



FIG. 1. (Color online) Schematic representation of a 5 QD 
device in a star configuration. The charging energies Ut and 
hopping matrix elements for each QD are indicated in the 
figure. Vl and Vr are the tunnel couplings to left and right 
electrodes, respectively. 

In the weak QD-electrodes coupling regime, the low en¬ 
ergy properties of the device are dominated by the ground 
state manifold of the isolated QD array which, as we show 
below, has spin S = (N — l)/2 and is ferromagnetically 
coupled to the electrodes. In the opposite regime, the 
magnetic moment of the central QD is screened through 
a regular Kondo effect leading to a local Fermi liquid 
of quasiparticles fIJ. The side-coupled QDs form a spin 
S = N/2 which is antiferromagnetically coupled to the 
central QD. To fully quench the associated magnetic mo¬ 
ment through Kondo screening, 2 S conduction electron 
channels would be needed. As we show below, there is 
a single conduction channel available and the magnetic 
moment is only partially screened. In both regimes and 
in the crossover between them, there is an asymptoti¬ 
cally free S = (N — l)/2 local magnetic moment at low 
energies. 

We calculate numerically the properties of devices with 
up to 6 QDs and obtain the expected logarithmic behav¬ 
ior in the thermodynamic and dynamic properties in the 
ferromagnetic and underscreened Kondo regimes and in 
the crossover between them. 

The rest of the paper is organized as follows: In Sec. 
[TT] we present the model for the N + 1 QD device. In 
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Sec. |III| we show that in the appropriate parameter 
regime, the low energy properties can be described by 


a S = (N—l)/2 ferromagnetic Kondo model. Section IV 


presents the numerical results, using the density matrix 
extension of Wilson’s numerical renormalization group 
(DM-NRG), for the spectral density of the central QD 
and the magnetic susceptibility. Finally, in Sec. [V] we 
present our conclusions. 


II. MODEL 


from the electrodes. We focus on the regime with a single 
electron on average on each QD and perform a Schrieffer- 
Wolff transformation (valid for: < 0, £( + U( > 0, and 
tn < min(|£^|, ee+Ui)) to decouple the high energy empty 
and double occupied states on each QD. 

Up to second order on the tunnel coupling term H t 
the Hamiltonian reads Hqd = Hj + Hyj + He- Here 
H\y describes density-density interactions between the 
central and the side-coupled QDs |21j and Hj describes 
the exchange interaction between the central QD and the 
side-coupled QDs: 


We consider N + 1 QDs, with a single relevant elec¬ 
tronic level on each QD [2t)j, coupled in a star configura¬ 
tion (see Fig. I]). The device is described by the following 
Hamiltonian W] 


H = H c + H t + H v + H el . (1) 


N 

Hj = ^ j(. S 0 ■ S(, 

1 =1 


( 6 ) 


where = YYaa' <4 ""f dio-' and er is a vector composed 
by the Pauli matrices er = (a x ,a y ,a z ). The antiferro¬ 
magnetic exchange couplings Je > 0 are given by 


Here He describes the electrostatic interaction on the 
QDs 


Ji = 8 tg 


U 0 + U e 


(Uo + U t f - 4(e 0 - 


(7) 


( ££ ^2 + Ueheyhei 

cr=t4 

where fii a = a/ ta di a is the electron number operator of 
the Uth QD, XJt > 0 is its charging energy, and a is 
the electron spin projection along the z axis. The level 
energy eg can be controlled experimentally via a gate 
voltage V g e, ei ~ —C g eV g e, where C g g is the capacitance 
of QD i with its corresponding gate electrode. 

N 

Ht = '^2"22 (ted\ rJ d 0a + h.c.'j , (3) 

(7 i—\ 

describes the tunneling coupling between the central QD 
(£ = 0) and the N side-coupled QDs. Finally, 

Hv = X! ( V ki' d l a C vk a + h.c)j , (4) 

V—L.R k,cr 

describes the coupling between the central QD and the 
left (L) and right ( R ) electrodes, which are modeled by 
two non-interacting Fermi gases: 


h c = J2 



Hel ~ ^ ^ e k c I/ ka Cl/ ka- (5) 

u,k,a 


III. LOW ENERGY HAMILTONIANS 


For the spin Hamiltonian of Eq. ([6]) two sublattices can 
be defined with no interaction between spins on the same 
sublattice, one formed by the central QD (sublattice A) 
and other formed by the side-coupled QDs (sublattice 
B). Under these conditions, the Lieb-Mattis theorem 
states that the ground state of the system has total spin 
S = — Sb\, where Sa ( Sb ) is the maximum possible 

total spin on sublattice A ( B ). This implies that the 
ground state of the isolated QD system will have total 
spin S = (N — l)/2. 

If we consider identical couplings Je = J between the 
QDs, the spin Hamiltonian reduces to an antiferromag¬ 
netic exchange coupling Hj —> J S 0 • Sq between the 
spin on the central QD and the total spin of the side- 
coupled QDs Sd = and can be readily solved. 

The eigenstates of Hj can be obtained adding the angu¬ 
lar momenta of the two sublattices L = So + Sd- Using 
JS 0 • S D = J(L 2 — Sp — Sq)/2 we obtain the correspond¬ 
ing eigenvalues: 

E L=s D ±k = {t(Sd + 1 ) + 5 )- ( 8 ) 

The minimum energy is obtained for the maximum possi¬ 
ble Sd = N/2 which results in L = (N— 1)/2 as expected 
from Lieb-Mattis theorem. The corresponding eigenstate 
is 

\L, L m ) = —j=== ^ \JSd + L m + d \i)\S D ,L m + \) 
- \JSd — L m + \ |t) | SoRm — 5 )^ , (9) 


A. Decoupled impurity 

We first analyze the nature of the eigenstates of Hqd = 
He + H t which describes the system of QDs decoupled 


where Sd = N/2, \a) and \Sd,S z ) are the eigenstates of 
Sq and Sd, respectively, with projections a and S z along 
the z axis, and L m is the projection of the total angular 
momentum along the same axis. 
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B. Ferromagnetic Kondo Hamiltonian 

In this Section we include the coupling to the electrodes 
and show that, in the appropriate parameter regime, the 
ferromagnetic Kondo model for an impurity spin S = 
(TV — l)/2 describes the low energy spin dynamics of the 
TV + 1 QD device. To that aim it is convenient to rewrite 
the coupling between the central QD and the electrodes 
[see Eq. @] as 

Hy = V (dl a c la + h.c)j , (10) 


terms for simplicity, and we have assumed N Je -C Uq. 
In the case that Je = J Jy, we expect the low en¬ 
ergy physics to be dominated by the ground state of the 
isolated QD array and we project Hi into its manifold: 
Hi —>• where | m) = \L,L m ). After the pro¬ 

jection, the first term of Hi turns into a constant that 
we discard. The projection of the second term can be 
evaluated using the Wigner-Eckart theorem based result 

m 

, IB , a M s o • Mm) , 

(m|So|m)= L ( L + Vj (m|L|m), (14) 


where V = ijYkv 1^1 2 and c ltr = y Yku v kvC k „a de¬ 
stroys an electron on the state T,) of the electrodes to 
which the central QD is coupled. 

To calculate the ferromagnetic coupling, we perform a 
second-order perturbation theory in Hy (Schrieffer-Wolff 
transformation) neglecting potential scattering terms 
and retaining only the magnetic coupling between the 
band electrons and the ground state manifold of the iso¬ 
lated QD array [see Eq. ([9])] . Following the notation and 
procedure of Ref. [T9] for Je = J we obtain an exchange 
coupling which is given by: 

J _ (L, L m + l|d 0 ! H Q n-E g — ^04. H Q n-E g ^0T^’ ^m) 

\j{L + ^) 2 — + I ) 2 

where L = S D - 1/2 = (N - l)/2, E g = E l = Sd -i /2 and 
we have neglected the energy of the conduction electrons, 
compared to the charging energies, since we are focussing 
our analysis on the low energy properties of the system. 
We finally obtain a ferromagnetic coupling 


which for the states of Eq. ([ 9 ]) is: 

1 


|S 0 |to') = - 


N + l 


ILIjn 7 ). 


(15) 


The minus sign of the proportionality coefficient be¬ 
tween So and L leads to a ferromagnetic coupling Jk = 
— Jy / (TV + 1) between s and L out of an antiferromag¬ 
netic coupling between s and So, and to recover Hj< as 
written in Eq. (12). Note however that to calculate Jy 


we have assumed IV J <C Uq, if this is not the case we ob¬ 
tain Jy = (jv+ 2 6 )j+ 2 £/ 0 an d recover the Kondo coupling 
of Eq. (11) for the case Je = J. 


In the general situation where the Je are not all iden¬ 
tical, the Hamiltonian Hj no longer commutes with Sp. 
The ground state of the isolated QD array [which is de¬ 
scribed by Eq. |6])] is, for a given L m , a superposition 
of the N states that have total spin angular momentum 
L = (N — 1)/2, of which a single one has Sd = N/2 while 
the remaining N — 1 states have So = N/2 — 1. It can 
therefore be written as: 


16K 2 

K ~ ~ {N + 1)(J(N + 2) + 2{7o) ’ 
for the Kondo Hamiltonian 


| L)gs = cosd|L, Sd = jf) +sin0| L,S D = jf - 1) (16) 

where 6 is a real number. Projecting Hi into the ground 
state manifold we obtain a Kondo coupling: 


H k = J K ~L-s + H el , (12) 

where s = Yaa' c -L“if“ c i!x' is the spin operator of the 
local state of the electrodes. This result, which extends 
the N = 2 analysis of Baruselli el al. to arbitrary TV, 
assumes that the ground state manifold of the isolated 
QD is the only relevant at low energies. As we will see 
below this is not the case if the coupling to the electrodes 
is strong enough. 

To show this, we perform the Schrieffer-Wolff transfor¬ 
mation taking Hy as a perturbation without projecting 
to the ground state manifold of the isolated QD array. 
We obtain the following Hamiltonian: 

N 

Hi=J2JeSo-S e + J v S 0 -s + H el , (13) 

e—1 

where Jy = 8V 2 /Uq is an antiferromagnetic coupling. To 
arrive to Eq. |l3|) we have discarded potential scattering 


( cos 2 9 sin 2 9 \ 

\ N+ 1 + N - 1 ) V ’ 


(17) 


which, as we show in the following section is always ferro¬ 
magnetic and attains its maximum absolute value for uni¬ 
form couplings (Je = J). The minimum absolute value of 
Jk is zero and is obtained in the limit where one of the Je 
couplings dominates over the others (e.g. Je^i/Ji —t 0). 


C. Proof of Jk < 0 


Here we show that the Kondo coupling of Eq. © 
is negative for arbitrary values of N and Je > 0. To 
calculate Jk from Eq. © only the coefficient 9 from 
the ground state wave-function of the isolated QD array 
[see Eq. (16)] needs to be obtained. It can be calculated 
diagonalizing the TV x N Hamiltonian matrix associated 
with the ground state subspace, which we managed to 
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do analytically for arbitrary N only in a few highly sym¬ 
metric cases. To show that the coupling is ferromagnetic, 
however, it suffices to obtain a lower bound on the value 
of cos 2 6. To that aim we propose a ground state wave- 
function of the form: 


Due to the positiveness of the we have 

= Yl a i +YU2 aia i - Yl a i’ ( 26 ) 

i / i i .1 /i i 


w = ( 18 ) 

i 

where the are real numbers and the 

= (19) 


which using the normalization condition of Eq. 
to and 


(22) leads 



(27) 


are “singlets” between the central QD and the i-th 
QD, with the rest of the spins parallel. Here, the ket 
|aoai • • • a^v) is a state in the Fock basis where the spin 
projections along the z-axis for the spin-1/2 on each QD 
are given by the ap. The |«S/) form a non-orthogonal 
basis of the subspace with total angular momentum 
L = {N — l)/2 where the ground state is to be found ac¬ 
cording to Lieb-Mattis theorem. This makes a variational 
calculation of the exact for the ground state wave- 
function. Applying Hj to the proposed wave-function 
the following linear combination of the | Si) is obtained 


Combining Eq. (27) and Eq. (251 we finally obtain 


Jk < 0, showing that the coupling is ferromagnetic. 
When one of the couplings dominates over the others 
—> 0 we have op = 1, o^i = 0 and Jk = 0. 
The maximum Jk in absolute value is obtained for sym¬ 
metric couplings Je = J for which at = 

Jk = — Jv/{N + 1). 


N(N+1) 


and 


D. Two-stage Kondo regime 


* l 


Jj T 


I Si), 


( 20 ) 

from which the variational energy E var = (\& | AT,/1H/) can 
be calculated 


E 


var 



ai<x j 

3 


1 

6 


yi a j 



that has to be minimized under the normalization con¬ 
straint: 


<w = \ 



( 22 ) 


Since all the coefficients in the quadratic form of Eq. 
© have the same negative sign, the a,; that minimize it 
will also have the same sign and can be chosen as positive 
(cti > 0) without loss of generality. The coefficient 6 can 
be calculated from the projection of T) on the state 


|L, = r,S D = JV/2) = /T^II^|s,), (23) 


which leads to 

cos 9 = (^>\L z =L,S d = N/2) = 


N + 1 


2N 2 > («) 

i 

Replacing the calculated cos 9 in Eq. © results in 

Jv 


Jk = 


N - 1 


HE 


Oii 


(25) 


The projection of Hi to the ground state manifold of 
the isolated QDs is justified if the energy gap between 
the ground state and the first excited state is much larger 
than the coupling energy to the electrodes. If we assume 
the opposite situation Jy J(, the low energy physics 
can no longer be described by the ferromagnetic Kondo 
Hamiltonian. For Jz —> 0, the model reduces to the usual 
spin 1 /2 Kondo Hamiltonian and the local magnetic mo¬ 
ment of the central QD is screened at temperatures below 
the Kondo temperature [21] 

Tk = ~j^ V PoJv exp[— 1/poJy], (28) 

where D is a high energy cutoff, and po is the local den¬ 
sity of states of the electrodes at the effective state IT!) 
evaluated at the Fermi energy. The local density of states 
at the central QD presents an Abrikosov-Suhl resonance 
of width ~ ksTft and the local low energy properties can 
be described by a Fermi liquid with quasiparticles having 
a renormalized mass oc 1 /T^ [T]. 

If we turn on a small Je = J <C the local spins 

on the side-coupled QDs couple antiferromagnetically to 
the local renormalized Fermi liquid at the central QD. 
This situation has been analyzed for Af = 1 using NRG 
and a slave-boson mean field approach in Ref. [25]. The 
main conclusion is that for J < hsT^ there is a second 
Kondo effect with a characteristic temperature: 

T* exp[—7rfc B T°/J], (29) 

which corresponds to a Kondo effect due to the antifer¬ 
romagnetic coupling J between a spin-1/2 and a Fermi 
liquid having an effective bandwidth D ~ k^T^- and a 
local density of states at the Fermi level po ~ 1 /nh^T^. 















A. Ferromagnetic to underscreened crossover 

(iV = 2) 
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FIG. 2. (Color online) Spectral density of the central QD of a 
device with two side-coupled QDs (N = 2). The parameters 
are Ue = U = 0.4, et = —U/2, A = 0.04-7T and te = t = 0, 
0.015, 0.02, 0.025, 0.03, 0.05, 0.07, 0.09, 0.11, and 0.13. Left 
inset: plot of Ao(w) -1 ' 2 as a function of ln(u;) to make clear 
the singular behavior of the spectral density at low energies 
for t A 0. Right inset: low energy detail of the spectral density 
for U=t = 0, 0.015, 0.02, 0.025, and 0.03. 


As we show below, for N > 1 the situation is anal¬ 
ogous to the N = 1 case in the J -C regime. The 
main difference being that for N = 1 the spin on the 
side-coupled QD is fully screened while for N > 1 the 
conduction electron channel available at the central QD 
is unable to fully screen the local magnetic moments at 
the side coupled QDs, leading to an underscreened Kondo 
effect. A numerical analysis shows that the characteris¬ 
tic temperature of the underscreened Kondo effect is well 
described by Eq. (29). 


Figure [2] presents the zero-temperature spectral den¬ 
sity of the central QD Ao(w,T = 0) in a N = 2 device, 
with identical QDs (Ug = U, el = —U/ 2), for different 
values of the interdot tunnel coupling ti = t. In this pa¬ 
rameter regime there is a single electron on average on 
each QD. For t = 0 we obtain an Abrikosov-Suhl reso¬ 
nance at the Fermi level associated with the screening of 
the magnetic moment of the central QD. The full width 
at half maximum of the resonance is ~ ksT^ and its 
height at the Fermi level Aq(uj = 0) = 2/7 tA according 
to Friedel’s sum rule [35] for the Anderson model. There 
are also two peaks at lo ~ Eq, and Uq+Eq associated with 
charge fluctuations. 

When a small coupling t between the QDs is included 
the spectral density develops a dip and vanishes at the 
Fermi level. The dip is associated with the underscreen¬ 
ing of the magnetic moment of the side-coupled QDs and 
leads to a complete suppression of the spectral density 
at the Fermi level at zero temperature. This dip can be 
interpreted in terms of a Kondo hole formed as a conse¬ 
quence of the strong coupling between the side-coupled 
QDs and the central QD. For N = 2 the vanishing of 
the spectral density at the Fermi level is contrary to the 
expectations from Fermi liquid theory and Luttinger’s 
theorem [5S] that predicts for this system in the wide- 
bandwitli limit |37j 


<> = o,r = o) = £ si, AjAl), (so) 

a 

where N a is the total occupancy per spin of the QD array 
[38] , In the absence of a magnetic field we have Nf = 
N\. and in the electron-hole symmetric situation we have 
N a = (N + l)/2, and 


IV. NUMERICAL RESULTS 


We calculate the thermodynamic and spectral proper¬ 
ties of the QD system using Wilson’s Numerical Renor¬ 
malization Group l2fiH321 with the density matrix exten¬ 
sion [33] and z-trick [34] to improve the accuracy of the 
low temperatures properties and the high energy resolu¬ 
tion of the spectral densities, respectively. The electron 
bath provided by the electrodes is characterized by the 
hybridization A(w) = Trpo(u>)V 2 where po(co) is the local 
density of states of the electrodes at the orbital | i/ji) to 
which the central QD is coupled. For simplicity we take 
A as energy independent and its support within the range 
[—D,D], where D is the half bandwidth of the conduc¬ 
tion electron band. In what folows we select energy units 
such that D = 1 and take the Fermi energy at ef = 0. 


Aq L (w = 0, T = 0) = 


2/ttA 


for odd N , . 

for even N ' ^ ' 


For t > 0 the spectral densities of Fig. [2] do not satisfy 
Eq. (31) indicating a non-Fermi liquid (NFL) behavior 
and a non-vanishing Luttinger integral [35]. In the left 
inset to Fig. [2]we plot [Ao(w)]~ 1,/2 as a function of ln(w) 
which shows a linear behavior at low energies consistent 
with a singular behavior of the form 


Aq(ui —> 0) ~ 


ln 2 (M /k B T o y 


(32) 


at low energies, as expected for the ferromagnetic and 
underscreened Kondo models [H HD] • The width of the 
dip increases with increasing t (see right inset to Fig. [2]) 
and suppresses the Abrikosov-Suhl resonance. For large 
enough t the Kondo effect associated with the screen¬ 
ing of the central QD does not develop and the system 
is in the ferromagnetic Kondo regime. The behavior of 
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FIG. 3. (Color online) Behavior of the low energy scale 
ksTo as a function of the interdot antiferromagnetic coupling 
J = At 2 /U for a system with two side-coupled QDs. Other 
parameters are as in Fig. [2] The lines represent the behav¬ 
ior of the Kondo temperature in the two correlated regimes. 
The dashed style line is the Kondo scale for a ferromagnetic 
model. The solid style line is a fit using the formula for the 
second stage Kondo temperature [see Eq. (291] with fitting 
parameters ci = 0.9, C 2 = 3.8. 


FIG. 4. (Color online) Magnetic moment squared of a triple 
QD device as a function of the temperature for different values 
of the interdot hopping t. All parameters are as in Fig. [2] 
except for the t = 0 curve which is not shown here. In the 
two stage Kondo regime {At 2 /U <S fcsT^) there is a two peak 
structure, while in the ferromagnetic Kondo regime {At 2 /U 3> 
k bTk) there is a single peak. The expected values of g 2 in 
four limiting cases (see text) are indicated on the right axis. In 
the whole range of values of t, the low temperature behavior 
is logarithmic as shown in the inset. 


the low energy scale ksTo as a function of J = 4 1 2 /U 
confirms this crossover from two-stage underscreened to 
ferromagnetic Kondo effects for J ~ ksT^. Fitting the 
numerical Aq{ui -a 0) using Eq. (32) we obtained To as 
a function of J which we plot in Fig. [3] The dependence 
of To on J for J <C ksT'x is consistent with an under¬ 
screened Kondo regime T 0 ~ T£ [see Eq. (29)]. For 


J ksTfr the temperature scale T 0 is consistent with 
a ferromagnetic Kondo regime with a coupling Jx given 
by Eq. dni. 


The underscreened and ferromagnetic Kondo regimes 
can also be identified through the behavior of the mag¬ 
netic susceptibility contribution \QD of the QD array as 
a function of the temperature. The magnetic susceptibil¬ 
ity XQD can be calculated assuming a coupling term: 


H B = H - B 


(33) 


between the magnetic moment of the QD array g = 
g e HB Sf and an external magnetic field B. For a free 
total spin L, the temperature dependence of the magnetic 
susceptibility follows a Curie law: 

C 

XQD = (34) 

where C = ( g e gB) 2 L{L + l)/3 is the Curie constant. In 
the general case, the spin in the QD is not free and sev¬ 
eral spin multiplets may be relevant at a given temper¬ 
ature. It is nevertheless useful to define a temperature 
dependent magnetic moment squared g 2 {T) = xqdT , to 
analyze the screening processes. The QD array contri¬ 
bution to the magnetic susceptibility is calculated, us¬ 
ing the NRG, from the fluctuations of the magnetization 


at zero field subtracting the contribution from the elec¬ 
trodes [MJ. In Fig. [i] we plot g 2 as a function of temper¬ 
ature. In the high temperature regime (fc^T U, t, A) 
all states of the QD array are equally probable which 
leads to g 2 ~ 3/8 {g e g B ) 2 , with each QD contributing 
~ 1/8 {g e gB) 2 to the magnetic moment. For J <C ksT 
g 2 {T) has two peaks as a function of the temperature. 
The high temperature peak at T ~ U/kB is associated 
with the formation of a spin-1/2 magnetic moment on 
each of the QDs which for three isolated QDs would lead 
to an increase in g 2 up to 3/4 {g e gs) 2 - The Kondo screen¬ 
ing of the magnetic moment of the central QD leads to 
a decrease in g 2 as the temperature is lowered which for 
J —> 0 would reach a plateau at 1/2 {g e gs) 2 associated 
with the two remaining uncoupled magnetic moments on 
the side-coupled QDs. The low temperature peak in g 2 is 
due to the coupling of the magnetic moments of the side- 
coupled QD (through the central QD) to form a triplet 
state with spin 5 = 1 [the magnetic moment squared as¬ 
sociated with a spin 1 is 2/3(geMs) 2 ]- For lower temper¬ 
atures the spin triplet is partially screened by the Kondo 
quasiparticles at the central QD leading to a residual 
magnetic moment g q = 1/4 {g e gs) 2 as T —► 0, due to an 
asymptotically free spin-1/2. 

As J increases {t increases) the two peaks in g 2 de¬ 
crease their amplitude and for J > Zc/jT/- the low tem¬ 
perature peak disappears. In the latter regime, the mag¬ 
netic moments of the QDs couple to form a spin 1/2 
at a temperature T ~ J/fcs and the two screening pro¬ 
cesses of the J <C ksTx regime do not take place. The 
spin 1/2, which is coupled ferromagnetically to the elec¬ 
trodes, decouples asymptotically as the temperature is 
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lowered m leading to a zero temperature magnetic mo¬ 
ment pi = 1/4 (g e PB) 2 - 

In the whole range of values of t ^ 0 explored the 
low temperature magnetic moment presents a singular 
behavior of the form 


i4 I i - 


i 


ln(T/T 0 ) 


(35) 


as expected from the Bethe ansatz solutions of the under¬ 
screened and ferromagnetic Kondo models [12H33]. This 
is illustrated in the inset to Fig. @ plotting (p 2 - nl) 1 
and a function of ln(T). Fitting the g 2 data with the 
expression of Eq. (351 we find pi ~ 1/4 {g e g B ) 2 and a 
temperature scale T 0 ~ T 0 . 


B. N dependence of the low-energy properties 

To reduce the computational cost of the numerical cal¬ 
culations we consider in what follows a weak interdot 
coupling regime where we can decouple the degrees of 
freedom associated with the empty and double occupied 
states in the side-coupled QDs and describe the coupling 
between QDs using Eq. ©>■ We also focus, for simplic¬ 
ity, our analysis on the electron-hole symmetric situation 
and Ji = J. We have checked numerically that our con¬ 
clusions hold for a wide range of parameters where these 
conditions are not met, breaking the electron-hole sym¬ 
metry and the symmetry of the exchange couplings. 

The analysis of the previous section can be extended 
for any number of side-coupled quantum dots. As it was 
shown in Sec. |III C'[ for a weak coupling of the central 
QD to the electrodes fc^T<C J, the low energy Hamil¬ 
tonian of the system is a ferromagnetic Kondo model for 
a spin S = ( N — l)/2. As we show below, in the oppo¬ 
site situation where J <C kgT^ this is not the case. As 
the temperature is lowered below T/k the magnetic mo¬ 
ment of the central QD is screened, and at lower tempera¬ 
tures (T < J /kg) a magnetic moment S = N/2 forms on 
the side-coupled QDs which is partially screened below a 
characteristic temperature ~ <C T^. 

Figure [5] presents the spectral density of the central QD 
for systems with up to 5 side-coupled QDs. The systems 
are in the two-stage Kondo regime and Aq(lu) presents 
the same qualitative features for all values of N. Two 
charge transfer peaks associated with the spectrum of the 
central QD are located at to = eo, U + eo and a central 
peak of width with a dip at the Fermi level of width 
~ ksTft. Due to the electron-hole symmetry considered, 
the spectral density vanishes exactly at the Fermi level 
which disagrees with the expectations from Fermi liquid 
theory for odd N. The ground state is, however, not a 
Fermi liquid for any N > 1. The right inset in Fig. [5] 
shows that the quadratic behavior expected for a Fermi 
liquid is only obtained for N = 1, for larger values of N 
the spectral density vanishes logarithmically as described 



FIG. 5. (Color online) Spectral density of the central QD 
Ao(w) for devices with N = 1,2,3,4, and 5 side-coupled QDs. 
The devices are in the two-stage Kondo regime J = 0.004. 
The other parameters are Uq = 0.4, eo = —0.2, and A = 
0.04-7T. Left inset: plot of Ao(cj ) -1 ' 2 as a function of ln(w) 
to make clear the singular behavior of the spectral density at 
low energies for N > 1. Right inset: low energy detail of the 
spectral density. 


by Eq. (321 which is confirmed plotting A 0 (uj) 1 / 2 as a 
function of ln(w) (see left inset in Fig. [ 5 ]). 

The spectral density in the ferromagnetic Kondo 
regime is presented in Fig. [6]for systems with up to 5 side 
coupled QDs. There is no Abrikosov-Suhl resonace in this 
case but a suppression of the density of states at low en¬ 
ergies. For w —» 0 a logarithmic behavior is obtained (see 
left inset in Fig. [6]) with a characteristic scale To that 
is shown in the right inset of the figure. Tq shows the 
expected behavior Tq ~ D/kg\/\poJK | exp(— 1/pqJk), 
with Jk given by Eq. ( fll| ), as a function of N which is 
presented as a dashed style line in the figure. 


The low energy thermodynamic properties show a dif¬ 
ferent behavior as a function of the temperature in the 
ferromagnetic and two-stage Kondo regimes. Figure [7] 
presents the magnetic moment squared p 2 (T) of the QD 
array as a function of the temperature. We also present, 
as a reference, the values of g 2 {T) for the QDs decou¬ 
pled from the electrodes (V = 0). As expected in both 
the ferromagnetic Kondo and the underscreened Kondo 
regimes, the low and high temperature limits coincide 
with the isolated QD array results. In the T —)■ 0 
limit both in the underscreened and in the ferromagnetic 
Kondo regimes the QD array has a magnetic moment 
pi = p 2 {T = 0) = S(S + l)/3 = ( N 2 - 1)/12 as ex¬ 
pected for a spin S = (N — l)/2 ferromagnetically cou¬ 
pled to the electrodes and to a spin S = N/2 under¬ 
screened by the conduction electrons. In the opposite 
limit kgT V, J(,Uq the QD array is effectively decou¬ 
pled from the electrodes and from each other and there is 
a contribution to the magnetic moment of N/4 from the 
side-coupled QD and of 1/8 from the central QD lead¬ 
ing to p 2 {T —► 00) = ( 2N + l)/8. In both regimes, p 2 
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FIG. 6. (Color online) Same as in Fig. [5] for devices with N = 
2, 3, 4, and 5 side-coupled QDs in the ferromagnetic Kondo 
regime, J = 0.4. Left inset: plot of Ao(oi) -1 ^ 2 as a function 
of ln(ta) to make clear the singular behavior of the spectral 
density at low energies for IV > 1. Right inset: N dependence 
of the low energy scale fcsTb- Numerical data (filled dots) and 
theory (dashed style line). 



FIG. 7. (Color online) Magnetic moment squared of the QD 
array for devices with N = 2,3, 4, and 5 QDs side-coupled to 
the central QD. a) Two-stage Kondo regime (J = 0.004). b) 
Ferromagnetic Kondo regime (J = 0.4). The other param¬ 
eters are as in Fig. [5] The dashed lines correspond to the 
magnetic moment squared for the isolated QD array (V=0). 


approaches the zero-tem pera ture limit in a logarithmic 
way as described by Eq. (35) with fig ~ ( N 2 — 1)/12 and 


Tg~T 0 . 

In the ferromagnetic Kondo regime fi 2 (T) deviates only 
slightly from the isolated QD array results. The coupling 
to the electrodes produces a small deviation which van¬ 
ishes logarithmically as the temperature is lowered and 
the QD array effectively decouples from the electrodes. 
In the underscreened Kondo regime, however, the be¬ 
havior of fi 2 (T) is more complex and the two Kondo 
screening processes can be observed as the temperature 


is lowered. For ksT < Uq a magnetic moment forms in 
the central QD, leading to an increase in fi 2 , which at 
lower temperatures is screened by the conduction elec¬ 
trons. For thermal energies lower than J, the spins on 
the side coupled QD form a magnetic moment associated 
to a spin S = N/2 leading to a peak in fi 2 that reaches 
~ N(N + 1)/12 for fcflT ~ J. The magnetic moment is 
partially screened leading to a logarithmic reduction of 
fi 2 as T —> 0. 


V. SUMMARY AND CONCLUSIONS 

We have analyzed the properties of a N + 1 QD dot 
device in a star configuration and showed that, as a func¬ 
tion of the tunnel coupling V between the central QD an 
the electrodes, there is a crossover between a ferromag¬ 
netic Kondo regime and a two-stage Kondo regime. In 
the weak coupling regime, the low energy properties are 
described by a ferromagnetic Kondo Hamiltonian for a 
spin S = (N — l)/2 impurity. In the strong coupling 
regime two Kondo effects take place successively as the 
temperature is lowered, a first stage were the spin-1/2 
of the central QD is screened by a regular Kondo effect 
and a second stage were a spin S = N/2 in the side- 
coupled QDs is partially screened and reduced to a spin 
S = (N — l)/2 in the ground state. In the full range 
of values of V / 0 the ground state of the system has 
spin S = (N — 1)/2 and the low energy dynamic and 
thermodynamic properties have a singular behavior. 

A numerical analysis of systems with up to 5 side- 
coupled QDs shows a singular behavior in the low en¬ 
ergy spectral density of the central QD which determines 
the linear transport properties of the device and is ac¬ 
cessible experimentally via spectroscopic measurements, 
in which one of the electrodes is weakly coupled to the 
central QD. We will present in a forthcoming publication 
the thermoelectric and magnetoelectric properties of the 
device that also show clear signatures of the singular low 
energy behavior of the device and allow to distinguish 
the different strongly correlated regimes. 

This QD devices would allow the observation of the so 
far elusive ferromagnetic Kondo effect for a spin S > 1/2 
and the underscreened Kondo effect for a spin S > 1. 
Devices with more than five QDs may be difficult to con¬ 
struct therefore limiting the maximum value of the spin 
in the impurity models. The models described could be 
also applied to magnetic molecules coupled to electrodes 
or arrays of magnetic atoms on metallic surfaces }45j . 
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